Vibrational detection and control of spin in mixed-valence molecular transistors 
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We investigate electron transport through a mixed-valence molecular complex in which an excess 
electron can tunnel between hetero-valent transition-metal ions, each having a fixed localized spin. 
We show that in this class of molecules the interplay of the spins and the vibrational breathing modes 
of the ionic ligand-shells allows the total molecular spin to be detected as well as controlled by non- 
equilibrium transport. Due to a spin-dependent pseudo John-Teller effect electronic transitions with 
different spin values can be distinguished by their vibronic conductance side-peaks, without using an 
external magnetic field. Conversely, we show that the spin state of the entire molecule can also be 
controlled via the non-equilibrium quantized molecular vibrations due to a novel vibration-induced 
spin-blockade. 

PACS numbers: 73.63.-b, 85.65.+h, 71.70.Ej, 85.75.-d 
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I. INTRODUCTION. 

In recent years, major experimental advances have 
been made in contacting and measuringsingle molecules 
in three-terminal transport junctions detecting the 
vibration d H IH , spin [4[ and magnetic properties [H, Q . 
Using the electrical gate to control the transport one can 
perform a substantial analysis of the transport processes, 
even when detailed microscopic information about the 
junction is lacking. Single-electron and coherent multi- 
electron tunneling processes (cotunncling, Kondo effect) 
governed by strong Coulomb and electron-vibration in- 
teractions allow detailed information to be extracted, 
as demonstrated recently for a single oligo-phenelyne 
molecule [1, Q. More detailed information can be ac- 
cessed when mechanical control of a gateable junction Q 
is possible. In the light of this progress of experimen- 
tal possibilities the interesting question arises how quan- 
tum mechanical states involving both electronic and me- 
chanical degrees of freedom may be detected and, per- 
haps, controlled in transport measurements. A particu- 
larly interesting aspect of single-molecule devices is the 
strong coupling of the electron current to the mechani- 
cal motion and the fully quantum mechanical character 
of this motion, resulting in Franck-Condon (FC) reso- 
nances [1, [Toj . In addition, in molecular systems with 
magnetic ions, the spin degree of freedom becomes im- 
portant [EH [0, El El El- Thus interestingmagnetic 
electro-mechanical effects are to be expected [15|, [l6[ . 

Mixed- valence molecules exhibit this interplay of quan- 
tum nanomechanics and spin-tronics. They are crucial 
also as a building block for supra-molecular devices and 
serve as benchmark for such systems. In a mixed- valence 
dimer, sketched in Fig.[TJ an excess electron can be local- 
ized on either of two equivalent metal ions with a local 
spin, thereby locally distorting the positions of the ligand 
atoms coordinating the ion. As the electron becomes de- 
localized over the molecule the distortion is "dragged" 
along coherently. As discussed in detail by Bersukcr 
and Borshch this results in vibronic mixing i.e. the 
molecular eigenstates do not allow an adiabatic Born- 



Oppenheimer (BO) separation of the nuclear and elec- 
tronic motions [l7|. The electronic and vibrational de- 
grees of freedom become entangled into vibronic states 
due to a spin- dependent pseudo Jahn- Teller (pJT) effect. 
In general, pJT mixing of electronic and vibrational de- 
grees of freedom cannot be neglected whenever electronic 
energy surfaces come close in energy (not necessarily de- 
generate as in the Jahn- Teller effect). The pronounced 
dependence of the derealization, and hence of the pJT 
mixing, on the total spin of the molecule arises due to 
local direct exchange interaction on the ions (Hund's 
rule), which favors the excess electron spin to be aligned 
with the ionic spin. Simultaneously, in such molecules 
this spin-dependent kinetic energy gain is responsible for 
the ferromagnetic double-exchange interaction [18;] which 
competes with other types of exchange interaction. 



In this paper, we present transport calculations for 
a model representative of a class of mixed-valence 
molecules. We demonstrate that a single-electron trans- 
port current can both detect and control the molecular 
spin due to the non-equilibrium nature of the vibrational 
motion. As indicated above, this does not rely on weak 
spin-orbit effects, but rather on strong direct, kinetic and 
double-exchange mechanisms. The pseudo Jahn- Teller 
dynamics shows up in pronounced sets of vibronic con- 
ductance peaks which depend on the spin-values of the 
molecular excitations. This provides a way to detect 
the spin without a magnetic field and probe the in-situ 
properties of a mixed-valence molecular transistor. Con- 
versely, we show that the electronic transport current 
induces non-equilibrium quantized molecular vibrations 
which drive a pronounced population inversion among 
the spin-states. Such a molecule can thereby be switched 
to a state with a well-defined charge and spin by ad- 
justing the applied voltages. This vibration-induced spin 
control arises from the interplay of spin and vibrations 
intrinsic to mixed-valence molecules, which may open up 
new possibilities for detection of mechanical motion and 
dissipation, magnetic switching and bistability. 
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FIG. 1: Sketch: Mixed-valence dimer trapped between two 
electrodes, L, R. By applying a bias voltage, Vt = fiL — 
(j,r, electrons can tunnel from one electrode to the other via 
the molecule. A third gate terminal (not shown), coupled 
capacitively to the molecule, effectively shifts the molecular 
energy levels. 



II. MODEL. 

We study the non-equilibrium transport properties of 
a well-established model for mixed-valence dimers from 
chemical physics, see [13, [3 and the references therein 
for detailed discussions. The model describes the sim- 
plest type of molecule exhibiting the spin-vibration in- 
terplay, consisting of two identical transition-metal ions, 
labeled by i — 1,2. We account for orbital-degenerate 
electronic states, \icr), where an excess electron with spin 
projection a is localized on ion i. Within the molecule, 
the electron can tunnel with amplitude t between the 
ions via a mechanically stiff bridging ligand. In addi- 
tion, the atomic groups that form a ligand-shell around 
each ion can vibrate along the local totally symmetric 
nuclear coordinate Qi (breathing mode). In the absence 
of an excess electron, this vibration around Qi = is 
assumed to be harmonic with frequency to. The excess 
electron, however, distorts the ligand-shell of the ion on 
which it resides, leading to a significant shift of the po- 
tential minimum of this mode (change in bond-lengths). 
In units of the zero-point motion energy of the vibra- 
tion, this shift is equal to \/2X, where A denotes the 
dimensionless electron-vibration coupling. In line with 
many experimental findings, we assume charging effects 
(Coulomb blockade) to be strong enough that only two 
molecular charge states need to be accounted for. The 
vibrational parts of the Hamiltonians H^ ih for N = 0, 1 
excess electrons on the molecule, written in the molecular 
vibrational coordinates Q± = (Q 2 ± Qi)j\[2, read 



i=± 



H° ih - XujQ + + \ujQ-(ni - h 2 ) 

+tY,( d L d 2a + h.C.) 7 



(1) 



(2) 



where hi = J2 a d\ a di a is the occupation operator of ion i 
and d\ a creates an electron in state \io). The symmetric 
coordinate, Q+, is the molecular breathing mode, which 
couples to the total excess charge, N, of the molecule, 
resulting in a shift by A of its potential surface along Q+. 



The resulting FC transport effects have been investigated 
theoretically [HI H3, El H3 and found experimentally 0, 
0, 0|- I n contrast, the anti-symmetric mode, Q-, couples 
to the internal charge imbalance hi — h 2 . Together with 
the intra-molecular tunneling, t, this results in the pseudo 
Jahn- Teller (pJT) effect [nj], i.e. the Hamiltonian © 
mixes electronic and vibrational states of the mode Q-. 

One of the crucial aspects in a mixed-valence molec- 
ular transistor is that the non-trivial dynamics of the 
Q- mode depends on the relative orientation of the local 
transition- metal ion spins. The molecular Hamiltonian 
for N — 0, 1 excess electrons reads 



- J(si • S 2 + Sr ■ s 2 ) - j hS 1 ■ Si. (3) 

i=l,2 

The intra- ionic Hund interaction, Jh, couples the spin 
of the excess electron, = \ ^Z a a i decree' di a i , to the 
spin of the transition- metal ions, Si, where er denotes 
the vector of Pauli-matrices and Sj = if no electron 
is present on ion i. Together with the intra-molecular 
tunneling, i, incorporated in Eq. ([2]), the intra- ionic Hund 
interaction results in a double-exchange splitting 2t$ of 
eigenstates with total molecular spin S fl8| ]: 



ts_ 
t 



S- 



251 + 1 



^-<1, 



(4) 



where Si = S 2 denotes the spin-length of the equiva- 
lent ions (Si,S 2 ). Result (H]) for the effective tunneling 
strength ts is obtained by expressing H 1 in the electronic 
basis of total spin eigenstates using vector coupling coeffi- 
cients. In the semi-classical limit of large ionic spins [l8j . 
Si ^> 1 (not considered further below), this reduces to 
ts/t = S/2S\ — cos(6>/2), where 9 is the angle between 
the two classical ionic spins. This makes clear that the 
kinetic energy which can be gained by the excess electron 
is maximal for parallel ionic spins due to the strong intra- 
ionic coupling and is suppressed by cos(6*/2), i.e. by the 
electron spin-eigenfunction component quantized in the 
direction of the ionic spin. Importantly, in mixed- valence 
molecules Jh is much larger than the other energy scales, 
Jh ^> |J|,w,t (typical values: Jh ~ 1 eV, J, t ~ 1-100 
meV, uj ~ tens of meV) [13, 0, • This scale separation 
derives from the intra-iomc origin of Jh (direct Hund 
exchange), in contrast to the intra-molecular processes 
involved in the exchange J and hopping t. We assume 
a ferromagnetic coupling, Jh > 0, i.e. a less than half 
filled ionic shell. The excitations where Si and Si are 
aligned anti-parallel can therefore be neglected. Eq. ([3]) 
also incorporates the intra-molecular coupling J of the 
spins of different ions. From hereon we take the length 
of the ionic spins to be Si = S 2 = 1/2. The Hamiltonian 
for the charged molecule consists of an S = 3/2 and an 
5 = 1/2 diagonal block, with 2S + 1 sub-blocks on the 
diagonal, given by 



Hi 



H; 



ib\t=t s 



±JS(S + 1) + const. 



(5) 
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This makes explicit the interesting property of mixed- 
valence molecules, that the strength of the pJT effect 
depends on the total molecular spin S due to the 
competition between the local distortion (coupling Xlo) 
and the effective derealization of the electron (energy 
ts). For N = 1 we numerically diagonalize the blocks 
separately for each spin state S to obtain the molecular 
eigenstates, which are not of the adiabatic BO form. For 
N = the states trivially factorize in BO form. 



assumed to be much smaller than the former, 7q <C 7o- 
The strength of the spin-allowed relaxation of course 
depends on the type of the molecular vibration mode and 
the junction substrate. For examples and discussion of 
very slow relaxation in the context of photon-tunneling 
in single- molecule junctions, see [25[ and the references 
therein. 



III. TRANSPORT. 

The interplay of vibrational and spin degrees of free- 
dom can be demonstrated for the basic transistor setup, 
sketched in Fig. [T] described by the transport Hamil- 
tonian H = Hmv + + H ICS . Here i?MV = 
Eat=o,i X n H n X n - aV g N, where X N projects onto 
states with charge N and a is the gate coupling. The 
electrodes r = L. R at electrochemical potential fj, r = 
jj,±Vb/2 and temperature T are described by 

H Ies = J2( €k - ^)cl ka c rklT , (6) 

r=L,R k.a 

i.e. Vf, is the bias voltage, and couple to the molecule by 

ht = Yl EE T '4<w +h.c. (?) 

r=L,Ri=l,2 k,a 

Here c^ fc(T creates an electron with spin a in state k of 
electrode r and TJ! is the amplitude for tunneling to ion 
i of the dimer. We assume a sequential arrangement i.e. 
T[ = T| = y/T/(2np) and T| = = 0, which may 
experimentally be favored by appropriate design of "clip- 
ping" ligands. Here T denotes the tunneling rate. 

In three-terminal molecular junctions, the tunneling 
rate can be as small as T « 0.1 meV and is often the 
smallest energy scale in the problem, enabling accurate 
spectroscopy of e.g. molecular vibrations [JJ. We focus 
on the regime of voltages and temperatures (T w 10r) 
where single-electron tunneling dominates the transport 
as is the case in many experiments 0, H, 13, S H, 0|- 
Using a kinetic (master) equation we calculate the 
non-equilibrium stationary state occupations of the 
molecule for each charge (N) and spin multiplet (5), 
keeping track of the symmetric vibrational (Q+) quan- 
tum number as well as the quantum number for the 
entangled state of the electrons and the anti-symmetric 
vibrations For weak tunnel coupling the transport 

rates can be evaluated in lowest non-vanishing order in 
Ht (i-e. Fermi's Golden Rule). Additionally, we account 
phenomenologically for relaxation due to coupling to 
a dissipative environment (e.g. substrate phonons). 
We assume an energy-dependent density of states (rate 
<y s ( v )(_E7) — . (E 2 /lu 2 )) for transitions between states 
with equal spin (7J) and different spin (7q). The latter 
relate to spin-orbit coupling effects and are therefore 




N=0 N=l 



FIG. 2: (color online) (a) Differential conductance (dl/dVt) 
for V g vs. V b (J = 2.9a;, t = 1.5w, A = 0.5, F = 9-W~ 4 lj, 7 <T = 
3.2 • 10~ 3 r, 7 g = 10~ 4 r, T = 10' 2 lu) (red: dl/dVt > 0, blue: 
dl/dVb < 0). The double-exchange coupling leads to a spin- 
dependent gap size of the vibronic spectrum (see arrows). The 
spectrum of S = 3/2 is harmonic (signalled by equidistant 
resonance lines of small energy separation), while the one of 
S = 1/2 is anharmonic (non-equidistant lines), (b) Sketch 
of energy spectrum: The spin-multiplets are split due to the 
intra-molecular coupling J. For N = 1, each spin-multiplet is 
split again by approximately twice the spin-dependent intra- 
molecular tunneling, ts, due to the pJT effect. Vibrational / 
vibronic excitations are omitted for clarity. 
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IV. RESULTS. 

A. Spin-dependent pseudo Jahn- Teller effect: 
Identifying spin values. 

The differential conductance as function of the ap- 
plied voltages is shown in Fig. Ela), using a set 
of parameters representative for mixed-valence dimers 
with weak electron-vibration coupling and ferromagnetic 
intra-molecular coupling: J = 2.9cj, t = 1.5a;, A = 0.5. 
The transport spectrum displays a number of sharp, well- 
separated resonance lines, which are "dressed" by many 
more lines with small separations. Two pronounced pairs 
of excitations appear, corresponding to the 5=1/2 and 
5 = 3/2 spin-multiplets which are split due to the intra- 
molecular coupling J (c.f. Eq. |5])). Each multiplet is 
split approximately by twice the spin-dependent intra- 
molecular tunneling ts, leading to a double-exchange 
gap which is reduced by a factor 2 for the 5 = 1/2 
state (ti/2 = £3/2/2, see also Fig. [2jb)). More gener- 
ally, the spin parameters can be determined from the 
ratios of the intra-molecular splittings giving ts/ts-i = 
(5 + l/2)/(5- 1/2). A central result of this work is that 
an independent check of this assignment of the spin is 
provided by the vibronic lines "dressing" the above exci- 
tations. This type of excitations arises when the shape of 
the vibrational potentials has a significant charge depen- 
dence [22J. Here their occurrence indicates a significant 
pJT mixing in the N — 1 charge state which changes 
the frequency and additionally induces anharmonicity in 
the effective adiabatic potentials. If the pJT effect is 
weak (as for 5 = 3/2), these potentials are approxi- 
mately harmonic in both charge states, and the spacing 
between these lines (that are due to transitions between 
excited vibrational and vibronic states) is even and equals 
the small frequency difference. For a stronger pJT ef- 
fect (5 = 1/2) the potential in the N = 1 charge state 
becomes anharmonic and the lines are unevenly spaced. 
Clearly, in Fig. Ufa) the "dressing" of the lower pair of 
lines is more evenly spaced than the upper set of lines, 
confirming the assignment of high spin state at low en- 
ergy. Here we merely note that a detailed analysis of the 
pJT transport resonances allows the electro-mechanical 
parameters ts,X,uj to be determined quantitatively, by 
reading off the voltage positions of the resonances. Also, 
we have invoked an adiabatic picture for the interpreta- 
tion, which has only an approximate validity, and some 
qualitative transport effects violate it. This is, however, 
not essential here, see (2(| for details. The spin identi- 
fication works very well for t > X 2 uj, i.e. when the pJT 
effect is weak to moderate (for | J| > t the resonances of 
the two spin-multiplets can be considered as separate). 
Thus the intra-molecular ferromagnetic coupling is re- 
vealed by the transport spectrum at zero magnetic field, 
by double-exchange and vibronic effects. Finally, we note 
that the energy average of the total 5 multiplets split by 
double-exchange follows J 5, providing a third, indepen- 
dent check of the spin value assignment. 



B. Current blockade and spin switching. 

For a wide range of parameters, the model exhibits a 
second, even more striking effect, exemplified for J = 
0.5a;, t = 5.0a;, A = 1.5 in Fig. (3Ja). At low energies, the 
current is stepwise reduced when going deeper into the 
sequential tunneling region leading to negative differen- 
tial conductance (blue lines in Fig.[3fa)). Simultaneously, 
the occupation of the molecular state with zero-spin and 
no vibrational excitations grows, reaching over 90% (see 
Fig.^b)). Within this region the current is strongly sup- 
pressed due to the pronounced population inversion that 
stabilizes the charge to N = and the spin to 5 = 0. This 
vibration-induced spin-blockade provides another indica- 
tion for the spin properties of the mixed- valence molecule 
and additionally allows the spin to be controlled electri- 
cally. The effect is readily understood by considering the 
non-equilibrium vibrations induced by the electric cur- 
rent. First we note that in the low-bias region where 
the spin-blockade occurs, the direct transition by elec- 
tron tunneling from 5=1— >5=l/2 is energetically 
not yet possible and the transition 5 = 0^ 5 = 3/2 is 
generally forbidden by spin selection rule (AS = ±1/2). 
Now consider an electron which has just enough energy 
to excite a vibrational, Q+, (or vibronic, Q-) quantum, 
when entering / leaving the molecule (N = <-> 1). If 
the molecule does not immediately relax it can accu- 
mulate more quanta in subsequent tunneling processes 
as sketched in Fig. G3c). Eventually, when a sufficient 
amount of vibrational energy has been accumulated a low 
energy electron can be assisted to excite the molecular 
spin-system. This tunneling process brings the molecule 
to a lower spin state with 5 = 1/2. Finally, the molecule 
can relax to the 5 = state by a single tunneling process 
in which the excess energy is dissipated into the elec- 
trodes. Now the molecule is trapped in a state with fixed 
charge N = and spin 5 = and the current is sup- 
pressed: neither the 5 = 3/2 states (due to the spin 
selection rule) nor the 5 = 1/2 state (due to the low bias 
voltage Vb < ( J + t) — 2V g and a strong Coulomb interac- 
tion on the molecule) are accessible. Fig. [2Jb) shows that 
the spin-blockade is lifted at higher bias Vb « ( J+t) — 2V g 
where the direct process back to 5 = 1/2 becomes possi- 
ble, thereby confirming the above mechanism. 

Clearly, the vibration induced spin-blockade is ex- 
pected to break down when the excited spin state is 
too high in energy to be reached or when escape pro- 
cesses from 5 = become dominant already at low bias 
voltages. The escape processes, however, have to change 
the spin by at least 1 quantum and are therefore para- 
metrically weak, since they relate to spin-orbit coupling 
(7o ^ 7o) or higher order tunneling processes. For in- 
stance, the phenomenological spin-flip relaxation which 
we included is responsible for the small remnant cur- 
rent in the spin-blockade region. Secondly, the block- 
ing state has to be reached: the single-electron trans- 
port current "pumps" the vibrational system (rate cx T) 
when the temperature is lowered below the vibrational 
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3: (color online) (a) Differential conductance (dl/dVb) 
for Vg vs. Vt (J = 0.5a;, t = 5.0a;, A = 1.5, r = 9 ■ 10 _4 a;, 
7o = 1.8 ■ 10 _2 r, 7 g = 10~ 4 r, T = 10' 2 uj). (b) Current 
vs. Vb (red solid) at V g = —0.2a; and probability of S — 
multiplet with no vibrational/ vibronic quanta excited (green 
dashed), (c) The energy spectrum of the N — 0,1 charge 
states and the spin-blockade mechanism. The line-style dis- 
tinguishes the total spin values S of the states. The longer 
lines denote a vibrational/vibronic groundstate whereas the 
shorter ones are excited by at least one such quanta. The 
molecule is "pumped" by a sequence of tunneling events, each 
time changing the charge and exciting vibrational and / or vi- 
bronic quanta until the S = 1/2 spin state is reached. From 
there, the molecule falls into the blocking state (5* = vi- 
brational groundstate) via another tunnel process. Since the 
transition to S = 3/2 is forbidden by spin-selection rules, 
S — cannot relax and the current strongly suppressed. 



frequency ui (preventing thermal relaxation) and when 
the tunnel coupling is sufficiently weak. It is thus cru- 
cial that the intra-molecular vibrations are not strongly 
damped. However, relaxation rates can compete with the 
transport rates without destroying the vibration-induced 
spin-blockade as long as the S = 1/2 state can still be ef- 
ficiently reached using the vibrations. For the results we 
present here this is indeed the case. Finally, cotunneling 
processes are expected to affect both the access to and 
the escape from the blocking state. A full calculation of 
this effect is possible in principle [27[, but is prohibited 
by the large number of states required here to describe 
the vibration-induced spin-blockade. Inspection of the 
numerically calculated rates to second order in T for a 
truncated spectrum indicates that the presented results 
are robust against perturbations due to higher-order tun- 
neling. Importantly, these processes can be suppressed 
by reducing the tunneling coupling by appropriate choice 
of connecting ligand groups. 

We note that in this paper we have discussed the case 
of ferromagnetic intra-molecular coupling J. However, 
the vibration-induced spin-blockade is generic and also 
occurs for anti-ferromagnetic coupling J < provided 
that t > | J| and in this case leads to a stabilized excess 
charge N — 1 and high spin S = 3/2. 



V. CONCLUSIONS. 

Using a representative model we have demonstrated 
that transport through a mixed-valence molecular tran- 
sistor entails an interplay of delocalized excess electrons, 
localized ionic spins, ligand-shell vibrations and Coulomb 
blockade. For this class of molecules transport-induced 
intra-molecular vibrations depend on the spin and their 
energy can be transferred to the spin-system at specific 
voltages, subsequently "locking" the spin, vibration and 
charge in a well-defined state. The generic model, an- 
alyzed here in a non-equilibrium situation, relates nat- 
urally to mixed-valence molecules [13] such as Ru 2+ ' 3+ 
complexes with pyridine organic ligands of the Robin- 
Day class II or III. The effects predicted here provide 
several bridges between nano-electromechanical systems 
(NEMS) and spintronic devices and applications in this 
direction can be envisaged. Clearly, the predicted spin- 
blockade effect will be sensitive to local magnetic fields, 
mechanical energy dissipation and spin-orbit effects and 
sensing applications involving these can be considered. 
Also, the blockade effect indicates slow transport dynam- 
ics: this may be used in the context of switching and 
bistable operation of molecular transistors by sweeping 
voltages non-adiabatically, with the new possibility of 
magnetic field control due to the involvement of spin. Fi- 
nally, from a chemistry perspective, transport measure- 
ments provide unique insight in the degree of electron 
derealization determining the key properties of mixed- 
valence molecules embedded in an electric circuit. Thus 
the investigation of complex mixed-valence systems as 
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devices 28J proves to be an interesting avenue in single- 
molecule electronics. 
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